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Abstract 

We discuss a remarkable property of an iterative algorithm for eigenvalue problems recently ad- 
vanced by Waxman that constitutes a clear advantage over other iterative procedures. In quantum 
mechanics, as well as in other fields, it is often necessary to deal with operators exhibiting both a 
continuum and a discrete spectrum. For this kind of operators, the problem of identifying spurious 
eigenpairs which appear in iterative algorithms like the Lanczos algorithm does not occur in the 
algorithm proposed by Waxman. 
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The Hamiltonian operator which describes a quantum mechanical system generally pos- 
sesses both a continuum as well as a discrete spectrum. A similar situation also occurs in 
other fields, such as theoretical population genetics jOj. In many cases one is only inter- 
ested in a few of the lower-lying bound states of the system. When only bound states are 
present iterative algorithms such as the Lanczos algorithm j2|_yield good approximations to 
the lower-lying eignstates with good convergence properties J 0, 0,!^ . On the other hand, 
the presence of the continuum leads to complications which can be circumvented but not 
without introducing spurious eigensolutions that need to be identified and ehminated 
Such spurious eigensolutions, however, do not occur in an algorihthm recently proposed by 
Waxman^l. We compare the two algorithms and demonstrate this behaviour in a simple 
numerical example. 

The presence of the continuum leads to complications in the Lanczos algorithm Qj. Find- 
ing a suitable start vector is by no means trivial ^ since the Lanczos algorithm can only be 
applied to states which are normalizable in the C? sense. For those operators which possess 
a continuum as well as a point spectrum, the space spanned by the bound state eigenfunc- 
tions is by itself certainly not complete and a suitable start vector should be composed only 
of components in the subspace spanned by the bound state eigenvectors. Usually the start 
vector is chosen from a complete set of analytic C? functions which define a space T . This 
space is in most cases not necessarily of the same dimension as the subspace spanned by 
the exact eigenvectors. On the other hand, if the Lanczos algorithm is applied with this 
choice for the start vector, the eigenpairs obtained will correspond to those of the operator H 
projected onto T . A subset of these eigenstates must correspond to the exact eigenpairs of 
the unprojected Hamiltonian operator since the exact eigenstates can be expanded in terms 
of the complete set of states which span T . The exact bound states can be identified and 
separated from the spurious bound states in the following manner j7|. After each iteration, 
for each of the converging eigenpairs (e//?, |e//3)), = |ef^— < eiplH^leiis > \ (where / is the 
iteration number) is calculated and a determination is made as to whether A is converging 
toward zero or not. For the exact bound states of H, A must be identically zero while the 
other spurious eigenstates states of the projected operator should converge to some non-zero 
positive value. Provided sufficient iterations are performed, it is possible, in this manner 
to identify uniquely the approximate eigenpairs which ultimately will converge to the exact 
bound states. 
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In the present note we wish to point out that this difficulty is avoided in a recently 
proposed iterative algorithm for determining the bound state eigenpairs of linear differential 
operators such as the Schrodinger Hamiltoniany]. This algorithm has many advantages 
not the least of which is its simplicity and an excellent convergence rate. The eigenpairs 
are determined as functions of the strength of the potential in the following manner. For 
simplicity consider a one-dimensional eigenvalue equation 

[-dl - XV{x)]u{x) = ~eu{x) (1) 

lim u{x) = (2) 

where = A > is the strength parameter of the attractive potential (AV(x) < 
and V(x)^ as \x\ oo) and the energy eigenvalue, — e (with e > 0), is negative and 
corresponds to a bound state. Using Green's method a solution to eq((T} is given by 



u{x) = A / G,{x - x')V{x')u{x')dx' (3) 
where the Green's function Ge{x) satisfies 

[-dl + e]G,{x) = 6{x) (4) 
lim G,{x) = 0. (5) 

Normalizing u(x) at an arbitrary Xref 

u(Xref) = 1 (6) 

allows A to be written as (see eq©) 

^ ^ J G,{x')v{x')u{x')dx' 

which can then be used to eliminate A from eqQ 

^ JZ.G,{x-x')V{x')u{x')dx' 
""^^^ J G,{x')v{x')u{x')dx' ■ 

Using equations (|7j) and (jH)), A can be determined as a function of e in the following manner. 
For a particular choice of e eq® can be iterated 

_ JZoG,{x-x')V{x')Ur.{x')dx' 
"""^'^^^ J G,{x')v{x')Un{x')dx' ^ ^ 



until it converges and A can then be determined from eq(I7)). Repeating for different values of 
e yields a set of different values of the potential strength A. When enough points have been 
determined, a simple interpolation procedure yields the dependence of e on A. Note that no 
diagonalization is required. In spite of the necessity of interpolating, the rapid convergence of 
the numerical solution of eqQ makes the present algorithm extremely viable. Furthermore, 
a proof of coverconvergence has been given and the algorithm can be extended for the 
calculation of excited states ^. 

In order to demonstrate that spurious solutions do not occur in the aforementioned algo- 
rithm we have performed the following simple calculation. An inverse Gaussian potential 

V{x) = e 2 

with half width of 2yJ {2ln{2)) has been constructed which does not support any excited 
bound states. The Lanczos algorithm has been used to determine the eigenstates of the 
corresponding Hamiltonian operator in one dimension using 

Ux) = {x\l) = i-y/'e-^' 

as the normalized start vector. After 18 iterations the Lanczos algorithm yielded the ground 
state at Cig i = -0.475917 plus a spurious state at cig 2 = 0.529612. In the case of the 
ground state Aig 1 = 0.0218906 while Aig 2 = 2.09673 clearly indicating that the excited 
state is spurious. The Waxman algorithm, using Ui{x) = 1 for the determination of the 
the interpolating function e(A), yielded the ground state at energy = -0.479203. Here the 
aforementioned iterations were repeated until A = 1 yielded a value to within 10"'^. When 
the Waxman algorithm was used to find the first excited state it did not yield a solution for 
A = 1. Hence no spurious solutions were obtained with the algorithm. Only for A > 1.35348 
did the Gaussian potential support at least one excited bound state . 

The spurious states arise in the case of the Lanczos algorithm because the resulting matrix 
representation of the Hamiltonian operator in the Lanczos basis corresponds to projecting 
it onto the space JF. The diagonalization of the resulting projected operator yields spuri- 
ous eigenpairs. In the Waxman algorithm, iterations in each step are performed with the 
Hamiltonian operator and no projection or diagonalization is required. Hence, ultimately 
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only the exact bound states are obtained and there are no problems with spurious states. 
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